RENEWAL CONVERGENCE RATES AND CORRELATION DECAY 
FOR HOMOGENEOUS PINNING MODELS 



GIAMBATTISTA GIACOMIN 

o . 

^ Abstract. A class of discrete renewal processes with super-exponentially decaying 

^S) ■ inter-arrival distributions coincides with the infinite volume limit of general homoge- 

' neous pinning models in their localized phase. Pinning models are statistical mechanics 

r- { , systems to which a lot of attention has been devoted both for their relevance for appli- 

cations and because they are solvable models exhibiting a non-trivial phase transition. 
The spatial decay of correlations in these systems is directly mapped to the speed of 
convergence to equilibrium for the associated renewal processes. We show that close to 
criticality, under general assumptions, the correlation decay rate, or the renewal conver- 
I gence rate, coincides with the inter-arrival decay rate. We also show that, in general, this 

. is false away from criticality. Under a stronger assumption on the inter-arrival distribu- 

tion we establish a local limit theorem, capturing thus the sharp asymptotic behavior of 
correlations. 
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O I 1- Introduction and main results 

. 1.1. Set— up and generalities. We consider the probability discrete density K{-) con- 

I centrated onN = {l,2,...}. We choose K[-) such that for some a > and some function 

■ L{-) which is slowly varying at infinity we have 

, n>N 

where we are using the notation ^ b]\f when limAr^oo otv/^a^ = 1- We assume that 
K{-) is aperiodic, that is gcd{n : K(n) > 0} = 1. We recall that a function L(-) defined 
on the positive semi-axis is slowly varying at infinity if it is positive, measurable and if 
limj^oo L{ct) / L[t) = 1 for every c > 0. We refer to [1| for the full theory of slowly varying 
functions, recalling simply that both L(t) and 1/L{t) are much smaller than (as t — > oo), 
and this for any 5 > 0. 

We point out that (jl.ip and aperiodicity are implied by 

Kin) ™ (1.2) 

Starting from K{-), we introduce a family of discrete probability densities indexed by 
6 > 0: 

Kb{n) := c(6)K(n)exp(-6n), (1.3) 
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and c{b) = 1/ K{n) exp{—bn) (of course c(0) = 1). 

Our attention focuses on the renewal process r(6) := {to(6), ri(6), r2(6), . . .} with inter- 
arrival law K{y{-), that is the process defined by ro(6) = and by the requirement that 
{rj+i(6) — Tj(6)}j=o,i,... is a sequence of IID random variables and P(ti(6) = n) = Ki,{n). 
Note that t(6) is an increasing sequence of almost surely finite numbers and it can be 
looked upon equivalently as a sequence of random variables (in fact, a random walk with 
positive increments) or as a random subset of N U {0}. With this second interpretation 
we introduce the so called mass renewal function, that is 

Ub{n) := P(nEr(6)), (1.4) 

so that Ub{n) is the probability that the site n is visited by the renewal. Note that Uh{0) = 1 
and, since K{-) is aperiodic, there exists ng > such that Ub{n) > for every n > uq. 

1.2. The Renewal Theorem and refinements to it. We now make an excursus in 
the general renewal theory on the integer numbers. We consider thus a general renewal 
process with tq = and with inter-arrival taking values in N. For this we introduce the 
notation F{n) := P(ri = n), while the mass renewal function is denoted by u{-). The 
classical Renewal Theorem (see e.g. [1]) says that, if F{-) is aperiodic, we have 

u{oo) := lim u{n) = G [0,1]. (1.5) 

n^oo E[riJ 

Much effort has been put into refining such a result. Refinements are of course a very 
natural question when E[ri] = -|-oo {e.g. [SI [10]), as well as if E[ri] < +oo. In the latter 
case sharp estimates on u{n) — u{oo) have been obtained for sub-exponential tail decay of 
the inter-arrival distribution, like for example in the case of F{-) = K{-) and K{-) as in 
(jl.2p (we refer to [13] and references therein). 

When instead the inter-arrival distribution decays super-exponentially, like for example 
if F{-) = Kb{-) with 6 > 0, general sharp results are harder to obtain. What can be proven 
in general in fact is that, if there exists ci > such that lim^^oo exp(cin)F(n) = 0, then 
there exists C2 > such that lim„^oo exp(c2n)|u(n) — u(oo)| = 0. However the precise 
decay, or even only the exponential asymptotic behavior (that is the supremum of the 
values of C2 for which the previous equality holds), in general does not depend only on 
the tail behavior of the inter-arrival probability. This is definitely a very classical problem 
[151 [Hj, and a number of results have been proven in specific instances (see e.g [3l I17j). 
We are now going to treat this point in some detail. 

1.3. On super-exponentially decaying inter-arrival laws. From the very definition 
of renewal process one directly derives the equivalent expressions 

n— 1 ^ 

= '^{Q}{n) +'y]u{j)F{n- j) and u{z) = ^— , (1.6) 

i=o 1 - ^(^) 

with the notation f{z) = Yl'^=o f i''^) (/(") is the z-transform of /(•)) and z is a complex 
number. Of course /(•) is a power series and \z\ a priori has to be chosen smaller than the 
radius of convergence, which, for the two series appearing in (jl.6p . is at least 1. 

As a matter of fact, we are interested (in particular) in the radius of convergence of the 
series 

oo 

A{z) := y (n(n) - ^(cx)))^" = --— -. (1.7) 

' ^0 E[ri](l-z) 
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If we assume that limsup„^oo exp(cn)F(n) < oo for some c > 0, the radius of convergence 
of F{-) is at least exp(c), however it is not at all clear that the radius of convergence of 
A(-) coincides with the radius of convergence of -F(-). In reality the problem does not 
come from the singularity at 2: = 1 (-F(l) = 1) since it is easily seen that it is removable. 
Notice also that, when F{-) is aperiodic, F{z) = 1 on the unit circle only if 2: = 1. However 
there may be other solutions z to F{z) = 1 for z within the radius of convergence of F{-). 
And it may even happen that A(-) can be analytically continued beyond the radius of 
convergence of F(-). Let us make this clear by giving two explicit examples: 

• F(l) = l-p, F{2) = p and F{n) = forn = 3, 4, . . . [p e (0, 1)). The radius of 
convergence of F{-) is 00, but A{z) = p/{{l +p){l + pz)) and therefore the radius 
of convergence of A(-) is 1/p, and in fact, by expanding A{z) around 2; = 0, we 
obtain u{n) — u{oo) = {—p)"'{p/ (1 + p)) for n = 1, 2, . . .. 

• F{n) = — p)/p, p £ (0,1). In this case the radius of convergences of F{-) 
is 1/p, but A{z) = p for every z, so the radius of convergence is 00 and in fact 
n(n) — n(oo) = for every n > 1. 

These examples show that the tail decay of n(-) — u{oo) may have little to do with the 
tail decay of the F{-): in particular, changing fine details of F{-) may have a drastic effect 
on the decay of n(-) — u{oo). For further examples of such a behavior see in particular [3], 
but also Section m below. 

The main purpose of this note is, however, to point out that, in a suitable class of 
renewal processes, motivated by statistical mechanics modeling (see Subsection II. 5p . the 
tail decay of u{-) — u{oo) is closely linked with the tail decay of the F{-). We are in 
fact going to show that F{-) = Kh{-), that is in the set-up of § 11.11 the decay rate of 
{ui,{n) — Ub{oo)}n is equal to the decay rate of Kf,{-), if b is sufficiently small. And under 
the stronger hypothesis (|1.2p we control the sharp asymptotic behavior of Uh{n) — Uf,{oo). 

1.4. Main result. With the set-up of § 11.11 we have the following: 

Theorem 1.1. Given K{-) call 6o(£ [0,oo]) the infimum of the values of b > such that 
there exists z satisfying 1 < |2| < exp(6) and Ki){z) = 1. 

(1) For every choice of K{ ) satisfying (jl.ip we have bo G (0, 00] and for every b £ 
(0, bo] we have 

limsup — log |nb(n) — nb(cxD)| = —b, (1-8) 

n^oo n 

while for b > bo we have 

limsup — log |nb(n) — nb(cx))| > —b. (1-9) 

n^oo ri 

(2) For every choice of K{-) satisfying ()1.2p we have that for every b G (0, 60) 

Ub{n)-Ub{oo) ~ (c(6) - 1)2 ' ^^-^^^ 

which implies 

lim — log {ub{n) — Uf,{oo)) = —b. (l-H) 
n— ►oo n 
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Remark 1.2. When there exists zq, 1 < \z()\ < exp(b), such that Kb(zo) = 1 (therefore 
b > bo) one can easily write down the sharp asymptotic behavior of {uj,{n) — Ufe(oo)}„ in 
terms of the values of zq with minimal \zo\. As a matter of fact one has 

limsup — log — Ufe(cxD)| = — logl^ol > —b, (1-12) 

n— >oo IT' 

but the sequence changes sign infinitely often and, in general, the superior limit cannot 
be replaced by a limit (see Section d] for details). In Section d] we also provide explicit 
examples showing that bo can be arbitrarily small by choosing K{-) suitably. In all the 
examples we have worked out the inequality in (jl.9p is strict (for every b > bo), but it is 
unclear to us whether or not this is a general phenomenon. 

The proof of Theorem ILlT l) can be found in Section [2] which is devoted to the study 
of Rb := l/limsup„ \ub{n) — Ub{oo)\^^^, which of course is the radius of convergence of 
Afe(-), and to establishing that bo is not zero. Theorem 11.1( 2) follows instead by a direct 
application of a well established technique : we detail this application in Section [3l We 
point out that the validity of the results in [7j go beyond the assumption (jl.2p . but we 
do make use of the regularly varying character of K{-) in establishing 6o > 0. A closer 
look at the proof of 6o > however shows that when ^jj^nK(n) < oo (cf. (j2.14p ) the 
regular variation property is used only marginally and in fact Theorem 11.11 holds also for 
a number of sub-exponential (cf. [4J) distributions K{-) beyond our assumptions. For 
example Theorem 1 1 . 1 1 holds also for K{n) = L{n)n'^ exp{—n'^), with g G M and 7 S (0, 1). 



1.5. Homogeneous pinning models and decay of correlations. What motivated, 
and what even suggested the validity of the results in this note, is the behavior near crit- 
icality of homogeneous pinning models. As it as been pointed out in particular in [9], a 
large class of physical models boils down to a class of Gibbs measures that, in mathe- 
matical terms, are just obtained from discrete renewal processes modified by introducing 
an exponential weight, or Boltzmann factor, depending on Mn{t) := |r n (0, A^]|. More 
precisely if P is the law of r and the latter is the renewal sequence with inter-arrival 
distribution K{-), we consider the family of probability measures {Pjv^^jjygj^ defined by 

dip 1 

-^(r) = ^ exp mNir)) , (1.13) 

with Zn^p the normalization constant. Then one can show ([6].[1H Ch. 2]) that the weak 
limit Foo,/3 of {Pat./sItvgn exists for every /3 G M (to be precise, this statement holds for 
every f3 assuming (|1.2p . but it holds also assuming only (jl.ip if /3 > 0). The parameter 
P actually plays a crucial role. In fact if /? < then r, under Poo,/3) is a transient 
renewal and it contains therefore only a finite number of points (this is the so-called 
delocalized phase). If instead /3 > then r, again under Poo,/3; is a positive recurrent 
renewal with inter-arrival distribution given by Kb{-), with b = b{f3) unique real solution 
of ^„ i('(n) exp(— 6n) = exp(— /3) (this is the localized phase). Note that if /? \ 0, then 
b\0. 

We point also out that it is not difficult to see that b coincides with the limit as N tends 
to infinity of (log Z]\i^p)/N and it is hence the free energy of the system [Til Ch. 1]. In [9] 
and, more completely in [HI Ch. 2], one can find the analysis of b{f3) as /3 \ 0. 

As a consequence r(6), for 6 > 0, does describe the localized regime of an infinite volume 
statistical mechanics system: if b is small, the system is close to criticality. The correlation 



CRITICAL PINNING AND RENEWAL CONVERGENCE RATES 



5 



length is a key quantity in statistical mechanics, see e.g. [U]. Moreover it is expected to 
scale nicely with P (or, which is equivalent, with b) approaching criticality, typically as P 
to some (negative) power, possibly times logarithmic corrections. The correlation length 
may be defined by introducing first the correlation function: 

P (m G r(6), m + n£ r(6)) - P (m G t(6)) P (m + n G r(6)) 



c(n) := lim 



°o Y^P (m G r(6)) (1 - P (m G t(6))) P (m + n G t(6)) (1 - P (m + n G r(6))) 
^(nGr(6)) ^ 



E[ri(6)]-lV ' ' E[ri(6)], 

(1.14) 

where we have used the Renewal Theorem. Then the correlation length is just one over 
the decay rate ^(6) of c(-): ^(6) := — 1/ lim sup„_,oo "-"^ loS 1^(^)1 therefore 

^(6) = — l/limsupn^"*^ log |tife(n) — nb(oo)|, (1-15) 

n— >oo 

SO that Theorem 11.11 guaranties that 

m ~° p (1-16) 

which roughly can be rephrased by saying that the correlation length, close to criticality, 
scales like one over the free energy. 

On physical grounds (jl.l6p . or rather the weaker form log ^(6) ~ —log 6, is certainly 
expected [9] . A proof of (|1.16p has been given in [H] by coupling arguments for the case 
in which K{-) is given by the return times of a simple random walk (and the proof is 
given also for disordered models). The result actually holds as an equality for every b 
(like the case presented in § 14.11 below: we point out that for a = 1/2 the distribution 
K{-) treated in § 14.11 coincides with the distribution of the returns to zero of a simple 
random walk in the sense that K{n) is the probability that the first return to zero of 
a simple random walk happens at time 2n). In general coupling arguments yield sharp 
results on the rate when suitable monotonicity properties are present (see in particular 
|16j): the returns of a simple random walk are in this class. In absence of monotonicity 
properties coupling arguments usually yield only upper bounds on the speed of convergence 
(and hence lower bounds on the rate, see [1] and references therein): in [19] a coupling 
argument is given for disordered pinning models and it yields in our homogeneous set-up 
that limsup^^Q log.^(6)/log(6) < —1, under the stronger hypothesis ()1.2p . 

We conclude this introduction with two important remarks: 

Remark 1.3. Some of the papers we have referred to (in particular [UITT]) aim at explicit 
bounds that hold for every n, possibly at the expense of sharp asymptotic results. Also in 
our set-up the question of obtaining more quantitative estimates, particularly when 6 \ 0, 
is important and relevant for the applications. 

Remark 1.4. The class of pinning models we have considered contains the so called 
(1 + d)-dimensional pinning models. The name comes from the directed viewpoint on 
Markov chains: if one considers a Markov chain S with state space Z*^, the state space 
of the directed process {{n, Sn)}n is 'Z^^'^. The renewal structure in this case is simply 
given by the successive returns to G Z'^ by S or, equivalently, by the intersections of 
the directed process with the line {(n,0) G Z^^'^ : n = 0,1,2,...}. This viewpoint is 
important in order to understand the spectrum of applications of pinning models. We are 
not going to discuss this further here, and we refer to [111 I20j . but we do point out that 
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precise estimates catching the order of magnitude of the correlation length in a class of 
(ci+l)-dimensional pinning models, i.e. Gaussian effective surfaces in a ((i+l)-dimensional 
space pinned at an hyper-plane, have been obtained in [5]. 

2. The radius of convergence of 

In this section we work in the most general set-up, i.e. we assume (jl.ip . Recall the 
definition of bo from the statement of Theorem II. li 

Proposition 2.1. Rf, < exp(6) and, for every choice of K(-), bo > and therefore 
Rf, = exp{b) forb£ {0,bo]. 

Note that this result implies (jl.Sp and (II. 9p . 
Proof. We are going to show that Rf, < exp(6) by making use only of Kiy(exp{b)) < oo and 
of the fact that the radius of convergence of Kf,{-) is exp(6). 

Of course we may assume that A;,(-) is analytic in the centered ball of radius exp(6), 
since otherwise there is nothing to prove. Let us suppose that Af,(-) has an analytic 
extension to the open ball of radius R > exp(6). From (|1.7p we immediately derive an 
expression for Ki,{z) in terms of Aiy{z), for \z\ < exp(6), and this gives the meromorphic 
extension of Kij{-) to the centered ball of radius R. However we know that the radius of 
convergence of is exp(6) and that < Yln^i^) < oo if \z\ = exp(6). So the 

singularity of Kh{-) cannot be a pole and therefore Kh{-) does not have a meromorphic 
extension. This implies that Ab(-) cannot be analytically continued beyond the centered 
ball of radius exp(6). 

The question that we have to address in order to complete the proof of Proposition 12. 11 
that is proving 6o > 0, can be rephrased as: do there exist two sequences {bj}j, bj \ 
and {zj}j, 1 < \zj\ < exp(6j) such that Kb{zj) = 1 for every j'? Of course, if this is not 
the case, Ki){z) / 1 if log |2;|(> 0) is sufficiently small. 

We make some preliminary observations: first, we may assume Q{zj) > 0, since if 
Kb{z) = 1, we have Kb{z) = 1 too. Then let us remark that, by writing zj = rj exp(i0j), we 
can pass to the limit in the equation Kb^ (zj) = 1: by the Lebesgue Dominated Convergence 
Theorem we have that every limit point (1,^) of {{f'j,9j)}j satisfies 

^ir(n)exp(m6i) = 1, (2.1) 

n 

which gives = hy aperiodicity. This tells us that, for b small, singularities have 
necessarily positive real part and small imaginary part (in short, they are close to 1). 
Moreover, by monotonicity, we see that the imaginary part cannot be zero (and therefore 
we assume that it is positive, since solutions come in conjugate pairs). 

Let us now assume by contradiction that there exists a triplet of sequences 

{{bjh, (2.2) 
tending to zero, with the requirements that > 5j < bj, 9j > for every j and such that 
Kb{exp(bj — 6j) exp{i9j)) = 1 for every j. Of course the triplet corresponds to the poles 
of the associated A;,^. (•) function at Zj = exp {{bj — 6j) + iOj). We are going to show that 
such a triplet does not exist since we are able to extract subsequences such that 



Kb{Bxp{bj - 5j) exp{i0j)) ^ 1, 



(2.3) 
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for every j in the subsequence. 

Let us consider the auxihary sequence of non- negative numbers By choosing 

a subsequence we may assume that this sequence converges to a hmit point 7 G [0, 00]. 

We consider first the case of q G (0, 1). We distinguish the two cases 7 < 00 and 7 = 00. 
If 7 < 00 we have the asymptotic relation 

Y: K{n) exp(-5,n) sin(0,n) ^'^"^ ejL{l/e,) H ^M-l')Ms) ^^^^^ 

n JO * 

that fohows from a Riemann sum approximation and the uniform convergence property 
of slowly varying functions [4, § 1.5] if the sum is restricted to OjU G (e, The rest is 
then controlled for small n's (n < e/6j) by replacing sin(3;) with x and using summation 
by parts which tells us that 'Yl,n=i nK{n) is equal to "^^=0 ^i^) ~ NK{N) and the latter 
behaves for large values of N as N^~"L{N)/{1 — a) [H § 1.5]. For large n's the rest 
is controlled by using \exp{—5jn)sm{6jn)\ < 1. Overall the absolute value of the rest 
is bounded by c6j'L{l/6j){e^~°' + e") for some c > 0, with c not depending on e, for j 
sufficiently large (for example, 9j < e) and l\2Ah follows. 

Observe that the left-hand side of ()2.4p is asymptotically equivalent to the imaginary 
part of Ki,{ex.p{bj — 5j) exp{i9j)), apart for the multiplicative constant c(6j) = 1 + 0(1) G M. 
The integral can be explicitly computed and it is equal to 

(1 + 7^)"''V(1 - a)sin(aarctan(l/7)) , (2.5) 

which is positive for every 7 G [0, 00), therefore for j sufficiently large (12. 3p holds (the 
definition of r(-) is recalled in Section H]). 

If 7 = 00 instead we write 

Y K{n) exp(-(5jn) sin(6ljn) = Rf + R> , (2.6) 

n 

with the sum for n < e/6j and R^" is the rest (0 < e < 7r/2 is a fixed positive constant). 
Setting Ss ■= sin(e)/e we have 

Rf > SeOj riK{n)expi-6jn) s,r{l - a)L{l/6j) (j-)^?- (2-7) 

To obtain (12. 7p we have used summation by parts, namely the identity: 



Y^ nK{n) exp(— Jj-n) = K{n) exp(— 5j(n+l)) — (1 — exp(— (5j)) nK{n) exp(— 5jn). 

(2.8) 



n=l ?i=0 n=l 

On the other hand 



< exp(-(5,/e,)5) Y ' exp(-(5,/0,)e) ^^^{9^/8^ , (2. 

n>e/ei 



therefore 

o— 1 / n \ o— 2 



< cexp{-{6,/9,)e) < c' exp (-(<5,/0,)e) , (2.10) 
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where c, c' are positive constants (we have exphcitly used the fact that, for every ci > 1 
and every C2 > there exists C3 > such that L{x)/L{y) < ci{x/yY'^ whenever x/y > C3 
[H Th. 1.5.6]). Therefore — > as j — > cxd and for j sufficiently large we have 

^ K{n) exp(-(5jn) s\n{ejn) > -SeT{l - a)L{l/6j)^6f , (2.11) 

2 j 

n ■' 

and then also in this regime (j2.3p holds. 



The marginal case of a = 1 and nK{n) = +00 is treated as follows. 
If a G [0, 00) for the step analogous to (12. 4p we split the sum according to whether 
9jn < e or 9jn > e. Summing by parts we obtain 

N N-l N , . 

5^ nK{n) = K{n) - NK{N) E ^ =^ ^(^)' (^.12) 

n=l n=0 n=l 

where in the asymptotic limit we have used 01 Prop. 1.5.9a] that guarantees that L(-) 
is slowly varying and that lining L{n) / L{n) = +00. From this we directly obtain that 
the first term in the splitting, i.e. the sum over Oju < e, is bounded below by a positive 
constant, depending on e and 7 (this constant can be chosen bounded away from zero 
for 7 in any compact subset of [0,oo)) times 9jL{l/5j). The rest instead is bounded, in 
absolute value, by a constant (independent of 7) times 9jL{l/9j), for j sufficiently large 
(just use I sm{9jn) exp(— (5jn)| < 1). Using once again L{n) ^> L{n) for large n, we obtain 
that ^„ K{n) exp(— 7jn) sm{9jn) > for j sufficiently large. 

If instead 7 = +cxd we restart from (12.60 and, by proceeding like in (12.70 and ()2.9p . we 
obtain that for j sufficiently large 

J2K{n)eM-S,n)sm{9,n) > ^SsL{l/5,) (^f)5j - 2 exp (-(5,/0,)e) L{l/9j)9j/s, 

(2.13) 

which is positive for j sufficiently large and the case a = 1 and nK{n) = 00 is under 
control. 

Let us now consider the case of a > 1, together with the case a = 1 and nK{n) < 00 
and note that in the latter case L(-) vanishes at infinity. 

In these cases for every 7 G [0, 00] we use the splitting in (j2.6p and for j sufficiently 
large we have 

J2K{n)expi-6jn)smi9jn) > -s^^j ^ nir(n) - -L{l/9j)9fe-'^ , (2.14) 

2 oi 

n n 

and the right-hand side is positive (again, for j sufficiently large). This concludes the 
proof of Proposition 12.11 □ 

3. Sharp estimates 

Throughout this section K{-) satisfies ()1.2p . we assume 5 > and we set Vub{n) := 
Ub{n) — Uh{n — 1) for n = 0, 1, . . . (iife(— 1) := 0). We also introduce the discrete probability 
density /if, on N U {0} defined by 

IXh{n) := Kb{n)/mb, (3.1) 
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with rrib := and Ki){n) := l^j>„ -f^feO)- Let us observe that 

rribfibin) = Kb{n) V exp(-6i) — -Kb{n), (3.2) 

K{n) exp(o) — 1 

and that this directly imphes the properties 

EUf^^iMn-j) n-.o. Mn + l) n;C-exp(-6). (3.3) 

/Xfe(n) fib{n) 
We point out also that from (jl.6p we get 



Vn„(z) = (f>b iiTbiz)) , with (?:);,(z) := ^— , (3.4) 

rribZ 

at least for \z\ < 1, like for (jl.Th . Of course the domain of analyticity of </>;,(•) is C \ {0} 
and if we observe that, by direct computation, we have 

^^b[z) = — -, 3.5 

mb{l - z) 

one can then extend the validity of (j3.4p to all values of z satisfying \z\ < exp(6) and 
1^1 < inf{|C| > 1 : KbiO = 1}. 

Proof of Theorem \l.lY 2). Let us choose h < b^. We observe that the two properties in (j3.3p 
are the hypotheses (a) and (/3) of Xh Theorem 1]. Hypothesis (7) of the same theorem, 
that is that jXbiz) converges at its radius of convergence (exp(6)), is verified too. Since 
h < bo, {iTb{z) : \z\ < exp(6)} C C \ {0}, i.e. the range of the power series /2fe(-) is a subset 
of the analyticity domain of (/>{,(•). Therefore [TJ Theorem 1] yields 

Vub{n) (p'b {fTbiewib))) = -7^7— ^^4S:;2 ' (3-6) 

(/x;,(exp(6))) rub 



and by (j3.2p we have 



Vub{n) ~ — — -2 — K{n)exp{-bn). 3.7 

c 6) - ly 



We conclude by observing that this yields 



Ubin) = - Vubij) j-T^Kin) exp(-6n) = - J^, (3.^ 



and the proof is complete. □ 

4. Some examples and further considerations 

Recall that T(z) := exp(— t) dt for ^{z) > 0, that r(-) can be extended as a 

meromorphic function to C and that T{z + 1) = zT{z) for z ^ {0, —1, —2, . . .} (therefore 
r(n) = (n — 1)! for n G N). Much of the content of this section is based on the fact that 
for /3 E R \ {0, -1, -2, . . .} and |x| < 1 we have 

n=0 

This is just a matter of realizing that for n > 1 
d" 

^(1 - x)-f = /3(/3 + 1) . . . (/3 + n - 1)(1 - x)-^-", (4.2) 
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and the formula is the Taylor expansion in x = 0. 

Since sign(r(/?)) = (-1)^/^11 for /? < ^ N) the first terms of the series in (gj]) have 
alternating signs, but for n sufficiently large the sign stabilizes and, by Stirling's formula 

T{x) "^-"^ exp(-x)x^-(i/2)^^ (4^3) 

one readily sees that T{n — a)/n\ Therefore, with the help of (j4.ip we can 

build probability inter-arrival distributions with the type of decay we are interested in and 
for which the z-transform is explicit. 

Remark 4.1. It is not difficult to see that one can differentiate, say j times, the expression 
in (j4.ip generating thus sequences which decay like (logn)-'/n^+° and that, for sufficiently 
large n, do not change sign. This provides examples involving slowly varying functions. 

Since we are just developing examples and that generalizations are straightforward, we 
specialize to the case of — /3 = a G (0, 1). 

4.1. The basic example. In this section we study the case of 

Kin) := —^7 ~ — — -. (4.4) 

^ ' -r(-a)n! -r(-Q) ^ ^ 

Note that YlnLi K{n) = 1 follows from (j4.ip . with (3 = —a, as well as, with reference to 
(fL3|l . c{b) = 1/(1 - (1 - exp(-6))") and 

(l-(l-zexp(-6))") 
""^^'^ - (l-(l-exp(-6))")- ^'-'^ 
In defining z" for a non integer, we choose the cut line {z G M : z < 0}. With this choice 
(1 — z exp(— 6))°^, and therefore Kii{-), has a discontinuity on the line {z E M : z > exp(6)}. 

We observe that, for every 5 > 0, Kf,{z) = 1 for \z\ < exp(6) only if z = 1, therefore 
Theorem 1 1 . 1 1 holds with bo = oo. 



Remark 4.2. In the special case under consideration, but also in all the other cases con- 
sidered in this section, one can obtain and go beyond Theorem 1 1.1 1 bv direct computations. 
In fact if we set q{z) := (1 — z exp(—b))" we have for \q{z)\ < \q{l)\ 

l-K,{z) a{z)-q{l) q{l) j^^W^) ) ' ^ ' ^ 

Now we set 

and we note that {q{z)y = (1 — zexp(— 6))-'" and therefore once again (j4.ip provides the 
expansion for {q{z)y if ja ^ N and the n-th term in the power series (of {q{z)y) behaves, 
as ra ^ oo, like cexp(— n6)n~^~''", c 7^ 0. Note that if ja G N the arising expression is 
just a polynomial and hence does not contribute to the asymptotic behavior of the series 
expansion. 

Finally, the series expansion ^„r("^)(n)z" of i? 

mi,') can be controlled by observing that 
this function is analytic in the centered ball of radius exp(?j) and by using the formula 

r^'^Xn) = ^ / dz = !^PizM j'^ (exp(5 + ^e)) exp (-m^) d^, (4.8) 
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where the contour in the middle term is (say) \z\ = r, for r € (0, exp(/3)), and the last term 
is obtained by letting r /" exp(5), using the fact that -Rm(exp(6 + «6')) is bounded. In fact, 
from the explicit expression and by construction, one readily sees that Rm (exp(6 + iO)) is 
smooth except aX6 = 27r/c, /c G Z, where it is C'-^'""''^)'^-! . By using the fact that n-th Fourier 
coefficient of a function is o{n~^), we see that A"^\n) = exp(— 6n)o(l/nl-("^"'"^)"-l ). 

The chain of considerations we have just made leads to an explicit expansion to all 
orders for exp{bn){uh{n) — Ub{oo)) as a sum of terms of the form Cj^j2n~^^~°'^'^ , for suitable 
(explicit) real coefficients Cj^^^ (ji,i2 eN). 

4.2. Singularities and slower decay of correlations. From the basic example one 
can actually build a large number of exactly solvable cases that display the more general 
phenomenology hinted by Theorem ll.il in particular that, in general, 6o < oo. 
For example, fix m S N and define 

^ ( r(n — m — q)/ (— r(— q) (n — m)!) for n = m + 1, m + 2, . . . [A^) 
I for n = 1,2, ... ,m. 

Note that this is nothing but the previous choice of K{-) translated to the right of m steps. 
Therefore 

^ _ (l-(l-zexp(-6))") 

(l-(l-exp(-6))") • ^^-'"^ 
Once again the radius of convergence is exp(6), but this time, in general, it is no longer 
true that one cannot find a solution zq to K{y(zQ) = 1 in the annulus 1 < \zq\ < exp(6). 

Let us choose a = 1/2 and let us first look at the case of m = 1. One can directly verify 
that 

1 



zo = + i/8exp(6) 1 - - exp(-6) - 3 < -1, (4.11) 



solves Kb{zo) = 1, that it is the unique solution (except the trivial solution zq = 1), and 
|zo| < exp(6) for b > bo with 

bo := log ^3/2 + V2-\JV2 + 5/4^ = 0.248399... (4.12) 

So, if 6 > 6o, since zq is a (simple) pole singularity of Ab(-) we can write 

ZqKI{zq) (1 - {z/zq)) 

with /(•) a function which is analytic on the centered ball of radius exp(6). Therefore 

Ub{n) - Ub{co) = ^ Zq-" + e(n), (4.14) 
and limsup„_^oj3(l/n) log |e(n)| = —b. 

Remark 4.3. Note that zq = —1 — exp(— 6)/4 + 0(exp(— 26)) for b large, so that the rate 
of converge of Ub{n) — Uoo{n) becomes smaller and smaller as b becomes large. 

Going back to (|4.9|) , for m larger than 3 one can no longer explicitly find all the solutions 
z to Ki){z) = 1. However we have the following: 

Proposition 4.4. For every b > and a £ (0, 1) one can find m G N such that if K{-) is 
given by (j4.9p then there exists a solution zq to Kb{zQ) = 1 with 1 < \zq\ < exp(&). 
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Remark 4.5. In general, once the solutions to Kf){-) = 1 of minimal absolute value (in 
the annulus {z : 1 < \z\ < exp(6)) are known, it is straightforward to write the sharp 
asymptotic behavior of ui,{n) — Uh{oo). For example if zq is a complex solution, then 
also its conjugate is a solution. If these have minimal absolute value among the solutions 
and if they are simple solutions, for a suitable (and computable) real constants ci and C2 
(|ci| + \c2\ > 0) we have 

Uh{n) - Ub{oo) koT" (ci cos (narg (zq)) + C2 sin (narg (zq))) . (4.15) 
An analogous formula is easily written in the general case. 

Proof of Proposition \4-4\ In reality, we are going to do something rather cheap, but we 
are actually proving more than what is stated: we are going to show that for every fe > 
and every r G (0, exp(6)) we can find an m such that there are m zeros of Kij{-) — 1 in the 
annulus {z : 1 < |z| < r}. 

Given 6 > 0, since the only solution z to 1 — (1 — 2;exp(— 6))") = is z = 0, then for 
every r S (l,exp(6)) we have 

1 - (1 -rexp(-6 + i6l))") 
e 1 - (1 -exp(-6))«) 

Therefore (recall (|4.10p ) |i^;,(2)| > r"^Xr, if \z\ = r. Therefore for m sufficiently large we 
have |ir;,(2;)| > 1 for \z\ = r: let us fix such a couple (m, r). Rouche's Theorem (e.g. [21 
p. 153]) guarantees that if / and g are analytic in a simply connected domain containing 
the simple closed curve 7 and if \ f{z) — g{z)\ < \ f{z)\ for z G j, then / and g have the 
same number of zeros enclosed by 7. Let us apply Rouche's Theorem with f{z) := Kiy{z) 
and g{z) := 1 — Kb{z) and 7 := {z : \z\ = r}, so that \ f{z) — g{z)\ = 1 < \ f{z)\ for z G 7, 
by the choice of m. But i^fe(-) has precisely m + 1 zeros (they are all in 0) and therefore 
also 1 — Kf,{-) has m + 1 zeros enclosed by 7. Of course 1 — Kh{-) has a zero in 1 and all 
the other zeros have absolute value in (I,?'). □ 
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